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I consider two-dimensional large- quantum chromodynamics with scalar 
quarks with particular emphasis on its strong coupling regime which has 
not been studied so far. Techniques necessary to deal with the infinitely 
^ oscillatory bound state wave functions in the strong coupling regime are 

4^ developed. I derive an estimate for the ground state mass and show that (1) 

the lightest hadron in the theory is massless and (2) the ground state mass is 
continuous across the transition between the weak and the strong coupling. 
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Two-dimensional quantum chromodynamics with scalar quarks is defined 
by the Lagrangian density 



C = -\Fl^,Fr + {DMiD.c^r - ra^^\ (1) 



Here the scalar field belongs to the fundamental representation of the 
SU(A^) gauge group, g and m are positive constants having dimension of 
mass and the covariant derivative is -D^^i = d<^i + igA''^-(f>j. For simplicity 
and symmetry I consider the theory with one quark fiavor. The large- 
theory is obtained by taking the limit N ^ oo and — )• while keeping 
g'^N = constant. 

The theory with fermionic quarks was developed by 't Hooft [1] who 
obtained an integral equation for the quark-antiquark bound state. The 
corresponding equation in the scalar theory can be obtained by basically 
retracing his derivation. In the light cone gauge A_ = the quark-gluon 
vertex is g[pi + P2) while the quartic coupling g'^ (f>* (f> vanishes. In the 
large- limit the dressed quark propagator is 



2k^k- — m? — g'^N/jT — g'^\k-\/TrX -\- ze' 

where A is the infrared cutoff. This leads to the following equation for the 
bound state of the scalar quark-antiquark pair: 

M-'<f>{x) = I — V + / <t>{x) 

\\lx{\ - x)j 

dy, (2) 




[x) - (f){y) [x + y){2 -x-y) 
IT Jo {x-yY 4^x(l - x)y{l - y) 



where M is the mass of the quark-antiquark bound state, cf) is the bound 
state wave function, and x is the light cone momentum scaling variable. 
This equation was first derived by Bardeen et al. [2] and subsequently also 
studied by other authors [3,4]. 

In the above references attention was restricted to the weak coupling or, 
equivalently, heavy quark regime of the large- theory. In this Letter par- 
ticular attention will be paid to the strong coupling regime of the theory, 
although some weak coupling results will be derived as well. It turns out 
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that the theory exhibits markedly different behavior in its strong and weak 
coupling regimes. In particular, the bound state wave functions in the two 
regimes have entirely different properties. I will derive an approximate ana- 
lytic expression for the ground state mass as a function of the quark mass and 
show that, contrary to what was previously believed, the lightest meson in 
the theory is massless. I also show that the mass of the ground state changes 
continuously as we pass from one regime to the other thus indicating that 
the spectrum of the theory depends continuously on the coupling (or quark 
mass). 

First note that there is an equivalence between light quarks and strong 
coupling (and between heavy quarks and weak coupling) in the large- the- 
ory. The equivalence is based on the fact that the bound state equation 
(2) involves a dimensionless parameter, - — o. When - — o < 1 we have the 

weak coupling/heavy quark regime and when > 1 we have the strong 
coupling/light quark regime. In the limit of infinite quark mass the theory 
describes free scalar particles and can be exactly solved. The bound state 
equation reduces to 

M''(f>{x) = 



x{l — x) 

whose solutions are cf>{x) = 6[x — |) and M = 2m. In the limit — )• oo, 
we have massless or, equivalently, infinitely strongly coupled QCD. As I will 
show the theory can be almost exactly solved in this limit and the ground 
state is massless. 

The first step in tackling the bound state equation (2) is to identify the 
proper Hilbert spaces in both regimes. This is done by studying the boundary 
conditions that the wave function must satisfy such that the leading singu- 
larities on the right hand side of Eq. (2) cancel. In the weak coupling regime 
the coefficient of the kinetic term is positive and the singularity cancellation 
dictates that the wave function must vanish at the boundary points and 1 
as x^ and (1 — x)^; here < /3 < 1/2 depends on the mass and the coupling 
constant as follows: 

m = TrptanTrp. (3) 

TT TT 

This equation maps the entire weak coupling regime, < < 1, onto the 
entire available range of fi between and 1/2 and has no solution in the 
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strong coupling regime. That is, it has no rea/ solution. I will adopt here the 
rather unorthodox point of view of allowing imaginary exponents. Indeed, 
one can show that the imaginary exponent wave functions x^^ and (1 — xf^ 
lead to the required singularity cancellations when > 1. In this case the 
equation for /3 is 

^ m^ = - 7r/3tanh7r/3. (4) 



TT TT 

Trm^ 



This equation maps the entire strong coupling regime 1 < < oo onto 
the imaginary-/? axis between and 0.38187. . . Hence the wave functions 
of hadrons made up of massless quarks must behave at the boundaries as 

^.0.38187... _ ^y0.38187..._ 

I conclude that in the strong coupling regime the Hilbert space consists 
of iv^(0, 1) functions of the form 

(P{x) = A{x)x'f^{l - xy^, (5) 

with Ai^x) behaving near the boundary as x^ and (1 — Taking into account 
the symmetry of the equation, Ai^x) can be shown to be either of the form 

oo 

A{x) = Y,a^x\\-x)\ aoT^O, (6) 
fc=o 

or of the form 

oo 

A{x) = {\-2x)Y^ai,x^{\-xf, aoT^O*. (7) 
fc=o 

The physics of the imaginary exponent wave functions is worth comment- 
ing on. When the quarks are heavy and the coupling is weak, the probability 
of finding them near the boundary, i.e. at very large relative momenta, is 
vanishingly small. They deal with the singularities by staying away from 
the regions where the singularities are. Light and tightly bound quarks, on 
the contrary, have finite probability to be found near the boundary. The 
way they kill the singularities is by having wave functions which oscillate 
infinitely rapidly in the singular regions! 

The Hamiltonian matrix element between the states %l) and ^ is 

(^, E^) = im^ - ^) r^^p^dx + Viv/^4iM.. 
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, g'N f^f^ {r{x) - r{y)){m - Hv)) {x + y){2-x-y) ^^^^ 

The functions (5) have unusual properties and the evaluation of the in- 
tegrals is subtle. For example, the absolute value of A{x)x^^^{l — xf^ is 
independent of /3 and equal to even though both the real and the 

imaginary part, A[x) cos(/3 In a;(l — x)) and A[x) sin(/3 In a;(l — x)), are highly 
oscillatory near the boundaries. In general, in the product t/j*[x)(f>[x) there 
is no vestige of the highly oscillatory behavior of its factors: it goes smoothly 
to non-zero values at the boundary points. This, of course, leaves the singu- 
larities in the kinetic term out of control. It is less obvious, but equally true, 
that double integrals involving products of the wave functions at different 
points t/j*[x)(f>[y) also blow up for the same reason. Apparently, even though 
the boundary behavior of (f>{x) cancels the singularities in the integral equa- 
tion (2), it does not do the same in the Hamiltonian matrix element (r/j, H(f>). 
This is unphysical as we cannot allow the Hamiltonian matrix elements to 
have singularities which are not present in the eigenvalue equation (2). 

To remove the singularities I propose an alternative definition of the prod- 
uct t/j*[x)(f>[y) which will "remember" that (f> and t/j are highly oscillatory 
functions. The following prescription gives correct and physically meaning- 
ful results: replace the product t/j*[x)(f>[y) by 

{nx)m)rea = k^*i^)Hy) + (9) 
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where 



^(x) = ^^{x){x''{l - xY' + x-''{l - x)-'') 

= ^A{x){x'^^+'\l - xY^^^'^ + x'^^-'\l - xY^^-'^) (10) 

and let e — )• at the end of the calculation. Note that since the tilde regu- 
larization commutes with the complex conjugation, the Hamiltonian remains 
hermitian under the regularization. With this prescription the integrals are 
well defined and can be computed. 

In both cases (6) and (7) it is the constant term which leads to the singular 
integrands in Eq. (8). In other words, if we can establish that the matrix 
element {(/),H(/)) exists for the first order wave function (/){x) = x^^{l — xY^ , 
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the higher order terms will not pose any new problems. The first order 
wave function is also relevant for the question of positive definiteness of 
the Hamiltonian and will provide us with an estimate of the ground state 
mass. Here I present the first order calculation. Higher order corrections 
and calculation of excited states will be reported elsewhere. 

One may wonder whether the Hamiltonian is positive definite in the 
strong coupling regime since the kinetic term is proportional to the quan- 
tity — which is negative. However, the kinetic term integral turns out 
to be zero for the wave function (j){x) = x^^{l — xf^ for all /3, i.e. throughout 
the strong coupling regime! This is a direct consequence of the regulariza- 
tion method (9) and it eminently makes sense. Note that as we approach the 
strong-weak boundary from the weak side, the kinetic term decreases and 
reaches the value zero at the boundary because its coefficient, the renormal- 
ized mass squared, becomes zero. From then on the kinetic term stays zero 
throughout the strong coupling regime - but now because the integral is zero. 
Since the remaining terms in (8) are manifestly positive the Hamiltonian is 
positive definite. 

The double integrals are evaluated by first regularizing them with the 
tilde products (9), computing the integrals, expanding the result in e and 
then letting e — )• 0. The result can be expressed in closed form in terms of 
Euler's psi function: 



This equation is also an approximate formula (more precisely, an upper 
bound) for the mass squared of the ground state meson. Together with 
m^(/3) given in Eq. (4) it defines parametrically the mass of the ground state 
as a function of the quark mass. This function is shown in Eig. f , to the left 
of the point = 

When the quarks are massless, the upper bound (f f ) for the ground state 
mass squared is 0.08 (in units ^-^), leaving little doubt that the ground state 




TT 



X -/3 + Im 



|(i - z/3)(V'(-l + 2z/3) + V(l - 2z/3) - V(-^ + - " ^Z?)) 



(12) 
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is exactly massless and its wave function very nearly 



(^(x) = - (13) 

It is amusing to observe that the absolute value of this wave function is equal 
to the exact ground state wave function of massless quarks in the fermion 
theory, namely 1. 

In order to compare these results with the weak coupling theory, I derive 
now the weak coupling regime analogue of Eq. (11). The Hilbert space now 
consists of the functions of the form 

A{x)x^{l-xy, (14) 

with A[x) behaving near the boundary points as and (1 — a;)°. Let us do 
again the first order calculation of the Hamiltonian matrix element which 
will give us the ground state mass as a function of /3. The double integrals 
can be expressed in terms of F functions by using an integral formula given 
in Ref. [2] and another formula derived by employing the Wilf-Zeilberger 
recurrence-finding algorithm [5]. The resulting Hamiltonian matrix element 
is 



TT r(i + 2/3)2 



o nWr , 2/3 + 1 r(| + /3)^ , ^^^l+j3)^. 
vr/3tanvr/3^ + ^ r(l + 2/3)^ + ¥TW{m^) J ^ 

Together with m^{fi) in Eq. (3) this expression defines parametrically the 
mass of the ground state as a function of the quark mass in the weak coupling 
regime. In Eig. 1, this curve is to the right of the point = The fact 
that both formulas join continuously at the point = where the strong 
and weak coupling regimes meet is a nice confirmation that our procedure 
for dealing with divergent integrals in the strong coupling regime is correct. 
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FIGURE 1 Mass squared of the ground state meson vs. quark mass squared. Both axes 
are in units of g'^ N / 1: . The weak coupling regime is to the right of the point = 1, the 
strong coupling regime to the left of it. 

In conclusion, the present work develops large- scalar quantum cliro- 
modynamics for light and strongly bound quarks. The wave functions in- 
troduced here and the techniques developed in order to compute with them 
make it now possible to do calculations in the entire range of quark masses. 
They can also be used in the full theory involving both scalar and fermion 
quarks in the strong coupling regime. As I will report in a forthcoming 
publication, I obtain massless composite fermions in that theory. Footnote 

* This leaves out the functions which, while going to zero at the boundaries, 
do so slower than any power. An example is provided by (lna;(l — x))~^ ■ 
This possibility is the reason for writing and (1 — x)^ elsewhere in this 
article. While I cannot exclude this possibility, in the present work I shall 
discard it. 
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